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Abstract 

In quantum detection theory the optimum detection operators 
must commute; admitting simultaneous approximate measurement of 
noncommuting observables cannot yield a lower Bayes cost. The lower 
bounds on mean square errors of parameter estimates predicted by 
the quantum-mechanical Cramer-Rao inequality can also not be 
reduced by such means. 
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Quantum detection and estimation theory has been developed within the 

conventional framework of quantum mechanics, one of the principal tenets of 

which is that only observables associated with commuting operators can be 
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simultaneously measured on the same system. It has been suggested that 

this formulation is too restrictive, that noncommuting operators can be at 

least approximately measured on the same system, and that to include this 

possibility may permit more effective detection, as measured by a lower average 
[4 5] 

Bayes cost. * We wish to show that no such improvement can be expected. 

The simultaneous measurement of noncommuting observables has been 
treated by Gordon and Louisell.^ In order to approximately measure certain 
such observables on a quantum-mechanical system S, it is made to interact for 
a time with a second system A, termed the apparatus. It was shown that a 
suitably defined ideal measurement yielding approximate values of the noncom- 
muting observables can be based on the outcome of measurements of commuting 
observables on the apparatus A, or more generally on both S and A. What we 
must therefore do is apply quantum detection theory — with its restriction to 
commuting observables — to the combined system S+A. 

Suppose we are to decide among M hypotheses , H 2 ,*.., H^. Under 

hypothesis Hj the density operator for the combined system at time t is 
S+A 

Pj (t) in the Schrodinger picture. If at an earlier time to the density oper- 

• S+A /si I 7 ] 

ator is Pj (to), the two operators are related by 

Qj-A Q+A 4- 

pj (t) = u(t, t 0 ) Pj (t 0 ) u (t, t 0 ), (1) 


with 


r t 


U(t, t 0 ) = exp 


-1 


Hdt'/fi 


( 2 ) 


where H is the Hamiltonian operator for the combined system S+A and ft is 
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Planck's constant h/2ir. The operator U is unitary; that is, with U its 
Hermitian adjoint, UU + equals the identity operator 1. 

Let { n j } be a set of commuting projection operators forming an M-fold 
resolution of the identity, 


M 



(3) 


j=l 

On the combination S + A we are to measure these M projection operators at time 
t, and if the k-th yields the value 1, hypothesis is selected as true.^^ 

The average cost is then 


C 


M M 



i=l j=l 


C ij ".I, 


(4) 


where is the prior probability of hypothesis Hj and is the cost of 

choosing H-^ when Hj is true. Let (Ilj(t)} be the projection operators that 
minimize C when the system S + A is observed at time t; we call these optimum. 
Then by (1) the operators 


JIj(t 0 ) = U + (t, t 0 ) Jlj(t) U(t, t 0 ) (5) 

f 

will minimize C when S + A is observed at time t 0 . Because of the unitarity 
of U(t, t 0 ) , the set { n j ( 1 0 ) } also forms an M-fold resolution of the identity 
into commuting projection operators, and the IT j (to) are optimum at time to- 
Since the minimization is carried out over all possible M-fold resolutions of 
identity, the minimum Bayes cost C m £ n must be independent of the observation 
time t. 

Now let us roll time back to an epoch to before the system S has come 
into contact with the apparatus A. In the Schrodinger picture this amounts to 
applying the inverse unitary transformation U + (t, to) to the state vectors of 
the combined system S + A. Because S and A are independent at this time tQ, 
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the density operators p^ must now have the factored form 
P^ +A (t 0 ) = P. S (to) p A (t 0 ), j = 1, 


( 6 ) 


Furthermore, as the apparatus A before the interaction has no information 

A 

about which hypothesis is true, p ( t q ) in (6) must be independent of j. 
The Bayes cost is now 


C = Tr 


M M 

-‘WEE'J C ±J p : S<to) "i 

i=l j=l 


(7) 


Since S and A are completely uncoupled, and the state of A is independent 

of which hypothesis is true, there is nothing to be gained by observing A. 

S A A 

The optimum projection operators n i (t 0 ) factor as IK (to) 1 , where 1 is the 

S 

identity operator for the apparatus A, and the set { II j (to)) forms an M-fold 


resolution of the identity 1 for the system S, minimizing the Bayes cost 

( 8 ) 


M M 

5 s ■ c u Trip / <t » ) "d- 

i=l j=l 


A .A 


Since Tr[p 1 ] = 1, the minimum value of C is also the minimum value of C 


in (7) and equals the time-independent minimum Bayes cost C 


min* 


The decision 


among the M hypotheses made at time tQ is based entirely on the measurement 

of commuting observables on system S. 

Similar considerations apply to estimating the m parameters 0 = 

c 

(0 1 , 02 ,«.*j 0 m ) of the density operator p°(0) of a quantum-mechanical system 

S. A version of the Cramer-Rao inequality sets lower bounds to mean square 

[3] 


errors of unbiased estimates of 0^, 02 >..*, 0 m * 


Let Xj be an operator 


whose measurement on S yields an unbiased estimate 0j of the j-th parameter; 

A 

0j must be an eigenvalue of Xj . Although in order to be measured simultaneously 
on the same system the operators X^ must commute, the analysis leading to the 
lower bounds given in [3] does not require commutativity of the operators X ^ . 
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For the class of commuting operators yielding unbiased estimates of the para- 
meters 0 there will exist lower bounds on the mean square errors, and those will 
be greater than or at least equal to the bounds derived in [3] . 

Again including the possibility of measuring noncommuting observables 
cannot lead to lower bounds smaller than those in [3] . In order to measure 
such operators even approximately, a measuring apparatus A must be allowed to 
interact with the system S, and according to Gordon and Louisell T s treatment 
of the process, commuting operators will at the end be measured on the combined 
system S + A. in the Schrodinger picture the density operator p^*^(0, t) 
for S + A will have a time dependence similar to that in (1) . 

Referring to (7) of [3] we see that the symmetrized logarithmic deriva- 
tives (SLD) Lj (t) appropriate for determining the Cramer-Rao lower bounds when 
the measurements of Xj are made at time t are related to those appropriate for 
measurements made at tQ by 

Lj (t) = U(t, t 0 ) Lj(t 0 ) U + (t, t 0 ). (9) 

Then (13) of [3] shows that the matrix A that sets the lower bounds is indepen- 
dent of the time t of observation, again because of the unitarity of the 
operator U (t , t 0 ) . 

Once more we move back to an epoch to before the system and the apparatus 

S+A S A 

have interacted. The density operator p (0, t 0 ) factors as p (0, tp) P (t 0 ), 

A 

where the density operator p (to) of the apparatus A is independent of the 
estimanda 0. The SLD operators for calculating the lower bounds are now the 
solutions of the operator equations 

8p S (e, t 0 )/38 n = \ [p S (e, t 0 ) L n s + L n s p S (6, t 0 )], (10) 

A 

and they act only on system S, commuting with p and all other operators on the 
apparatus A. When taking the trace over the states of A to form the elements 
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of the matrix A, the density operator p is replaced by 1, and the lower 

S 

bounds depend only on p (0, to). Thus the apparatus A cannot help estimate the 
parameters 0 of S with smaller mean square errors than the lower bounds cal- 
culated by the quantum-mechanical Cramer-Rao inequality as applied to the 
density operator of system S alone. 

In [3, p. 238] lower bounds were calculated for unbiased estimates m x and 
m^ of the components of the complex amplitude y = + iny of a simple harmonic 

oscillator, which might represent a mode of the field in an ideal receiver in 
the presence of thermal noise. Those bounds are 

Var fi x > j (N + y) , Var m y > y (N + y) , 

where N is the mean number of noise photons. The noncommutativity of the SLD's 
L x and Ly used to derive these bounds does not invalidate them. It can be 
shown that if the mode is coupled with an ideal amplifier whose gain is high 
enough to raise the oscillator variables to the classical domain where they 
commute, error variances 

Var m x = Var my = y (N + 1) 

can be attained [8] . It is unknown whether commuting operators can be found 
whose measurement will yield unbiased estimates m x and my with variances lying 
between ~ (N + y) and y (N + 1) . 

The measurements we need to make on a quantum-mechanical system S for 
testing hypotheses or estimating parameters will always have to be effected by 
means of an auxiliary apparatus A, and this apparatus, subject to thermal and 
quantum fluctuations of its own, will ordinarily introduce additional random 
uncertainties. Each measurement procedure will have to be analyzed to determine 
what error costs it entails. Detection theory and estimation theory seek 
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lower bounds on these costs, and in doing so they minimize with respect to the 
entire class of possible detection or estimation operators that can be applied 
to the system. The resulting bounds are independent of the time of observation, 
and they cannot be reduced by using any auxiliary apparatus that initially 
posseses no information about the state of the system. 
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